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ABSTRACT 

When averaged over large scales, star formation in galaxies is observed to follow the 
empirical Kennicutt-Schmidt (KS) law for surface densities above a constant threshold. 
While the observed law involves surface densities, theoretical models and simulations 
generally work with volume density laws (i.e. Schmidt laws). We derive analytic re- 
lations between star formation laws expressed in terms of surface densities, volume 
densities, and pressures and we show how these relations depend on parameters such 
as the effective equation of state of the multiphase interstellar medium. Our analytic 
relations enable us to implement observed surface density laws into simulations. Be- 
cause the parameters of our prescription for star formation are observables, we are not 
free to tune them to match the observations. We test our theoretical framework using 
high-resolution simulations of isolated disc galaxies that assume an effective equation 
of state for the multiphase interstellar medium. We are able to reproduce the star 
formation threshold and both the slope and the normalisation of arbitrary input KS 
laws without tuning any parameters and with very little scatter, even for unstable 
galaxies and even if we use poor numerical resolution. Moreover, we can do so for 
arbitrary effective equations of state. Our prescription therefore enables simulations 
of galaxies to bypass our current inability to simulate the formation of stars. On the 
other hand, the fact that we can reproduce arbitrary input thresholds and KS laws, 
rather than just the particular ones picked out by nature, indicates that simulations 
that lack the physics and/or resolution to simulate the multiphase interstellar medium 
can only provide limited insight into the origin of the observed star formation laws. 

Key words: galaxies: evolution — galaxies: formation — galaxies: ISM — stars: 
formation 



1 INTRODUCTION 

The efficiency of star formation (SF) and its environmental 
dependence are fundamental, but poorly understood, ingre- 
dients for models of the formation and evolution of galaxies. 
Fortunately, for many purposes, the models do not require 
a detailed understanding of SF in order to make progress. 
Our lack of understanding can often be bypassed by using 
empirical scaling relations which give the efficiency of SF av- 
eraged over scales that are large compared to those of indi- 
vidual star clusters. Such scaling relations can be built into 
semi-analytic models and numerical simulations of galaxy 
formation, allowing the SF process to be treated as a black 
box calibrated to reproduce the observations. 



* E-mail: schaye@strw.leidenuniv.nl 
f E-mail: caius@strw.leidenuniv.nl 



Observationally, a combination of two SF laws is known 
to describe th e efficiency of SF reasonably well for nearby 
galaxies (e.g. iKennicuttl Il998al ): a SF threshold and a 
Kennicutt-Schmidt (KS) law: 
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where E* is the rate of SF per unit area and per unit time, 
Eg is the gas surface density, and E c is the threshold surface 
density for SF. It is import ant to distingui sh the empirical 
KS law, E* oc from the ISchmidtl (jl959h law, p* oc pg s , 
which is often assumed in theoretical models. We will reserve 
n exclusively for KS (i.e. surface density) laws and will use 
ns to denote the power-law index of Schmidt (i.e. volume 
density) laws. 

The existence of a surface density t hreshold f or SF 
has long been explained in terms of the iToomrd ([l964) 
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criterion for instability (e .g. iQuirkl [l972l : IKennicuttl Il989: 
iMartin &; Kennicuttl l2Q0lh , usually under the assumption 
of a constant velocit y dispe r sion. This explanation was, 
however, criticized by Sch avd ([20041 . hereafter S04; see also 
Schaye 2007), who argued that the 2-D Toomre analysis is 
not applicable to the warm interstellar medium (ISM) in 
the outer parts of disc galaxies and that it fares less well 
in comparison to observations than a simple constant sur- 
face density threshold, even allowing for a n ad-hoc re- s caling 
of the critical Toomre Q-parameter, as in IKennicuttl |l989) 
and subsequent studies. 

The fact that a constant threshold of E c ~ 10 Mq pc -2 
works well empir i cally, had been k nown for some time (e.g. 
lGuiderdonilll987l : ISkillmanl fl987h . [M megreen fe Parravanol 
( 19941 ) emphasized that the need for a cold gas phase im- 
poses a pressure threshold. S04 showed that the transi- 
tion from the warm to the cold gas phase triggers gravi- 
tational instability on a wide range of scales. He demon- 
strated that this thermo-gravitational instability can ac- 
count for the observed SF threshold, which he predicted 
to fall in the range E c ~ 3-10 Mq pc -2 (corresponding 
to a critical pressure P c /k ~ 10 2 -10 3 cm K and a criti- 
cal volume density nu ~ 10 _2 -10 _1 cm -3 ), with weak de- 
pendencies on the metallicity, UV radiation field, turbulent 
pressure, and the mass fraction in gas. Conversely, below 
the threshold the presence of ionizing radiation keeps the 
gas warm (T - 10 4 K) and stable. The model can there- 
fore also account for the observed constant velocity disper- 
sion of ^ 8 kms" 1 in the outer parts of disc galaxies (e.g. 
iMeurer et allll996l ). Several recent studies have found that 
the S04 criterion can account for high- r esolution observa- 
tions of SF in nearby disc galaxies (e.g. |de Blok fc Waited 
l200d : lAuld et all l2QQ6|) and of the format ion of star clus- 
ters in tidal arms" (|Mavbhate et al.ll2007h . The study by 
Ide Blok fc Walterl (l2006h is particularly relevant because it 
avoids azimuthal smoothing, which is inappropriate for lo- 
cal thresholds and which has often confused comparisons 
between theories of local SF thresholds and observation. 

For E g ^> E c , many studies have found that the simple 
KS law, with n ~ 1-2, can well represent the global, disc- 
averaged star formation rat e (SFR) over at least 5 decades 
in E g (see IKennicuttl 1 99 8al for a review) . IKennicuttl (|l998bh 
founcQ 

(1.4±0.15) 



E. = (2.5±O.7)xlO- 4 M yr- 1 kpc- 2 ( r ^^n 

as the best fit for his sample of 97 galaxies spanning 5 orders 
of magnitude in E g , ranging from quiescently star- forming 
disc galaxies to starbursts. The KS law can also be expressed 
in terms of the gas consumption time-scale 



t g EE g^A-^lMQpC- 2 )^- 



(4 ± 1) x 10 9 yr 



(3) 



1 Mq pc 



-(0.4±0.15) 

• (4) 



1 The normalization assumes a Salpeter initial mass function over 
0.1-100 M q. Note that the e rror bars on the normalisation and 
slope of the IKennicuttl |l998b) result must be strongly correlated, 
since all his data points are on one side of the pivot point of 
1 Mq pc - 2 . Changing one parameter by itself, will therefore lead 
to results that are inconsistent with the observations. 



The same law is also found to be a good description 
of the relation between the SFR and gas surface den- 
sity when they are azimuthally av e raged in radial bins 
(e.g. IZhang. Fall, fe Whitmord l200ll ; IWong fe Blitzl 120021: 
Komugi et all 120051 : feoissier et al.1 120061 : ISchuster et~aD 
20071 ) and even when they are measured in local cylindrical 



(2) 



bins, althoug h in that case the coeffi cient A may be some- 
what smaller ([Kennicutt et al.ll2QQ7h . 

Note that the agreement between the azimuthally av- 
eraged and globally averaged (i.e., averaged over an entire 
galaxy) SF laws is surprising. Since the observed KS law is 
non-linear (i.e. n 7^ 1), the global law should only equal the 
local one if individual galaxies have constant surface densi- 
ties. The fact that the same law works in both cases therefore 
implies that, in terms of their SF properties, galaxies have 
effectively constant surface densities. This will be approxi- 
mately true if the global SFR is dominated by the nucleus. 

The KS law depends on the total gas surface density, 
but several studies have claimed that the exponent n ^ 1.4 
is only valid for molecular gas and t hat including atomic 
gas results in lar ger exponents (e.g. IWong &~B litz 2002; 
iHever et al.ll2004h . However, the HI surface density is ob- 
served to saturate at Em ~ 10 Mq pc -2 , which can be ex- 
plained as the result of the formation of a cold gas phase 
in clo uds with surf ace densities that exceed this critical 
value (jSchavdl2001bh . Such clouds will become unstable and 
will eventually become mostly molecular, unless the process 
is halted by SF and associated feedback processes. Since, 
not coincident ally (S04), this surface density is similar to 
the threshold for SF, H I will contribute significantly to (or 
even dominate) the gas surface density for densities below or 
around the threshold. It is therefore not surprising that stud- 
ies focusing on this low density regime, where a significant 
fraction of the gas at a fixed radius will have a surface den- 
sity below the SF threshold, find steeper KS laws. Empirical 
prescriptions that take into acc ount the molecular fractio n 
may work better in this regime (Bli tz &; Rosolowsky||2006h . 
This does, however, not undermine the validity of the KS 
law for Eg ^> E c . 

Many theories have been proposed to explain the ob- 
served KS law (e.g. |Silkl|l997l: iTanl l2000l : lElmegreenl 120021 : 
iKrumholz feM cKee 200% iTutuk ov 2006). Most rely in some 
way or another on the assumption that the SFR is propor- 
tional to the gas density over the dynamical time-scale, 



p* = e 



tdyn 



■ eG 1/2 pl /2 , 



(5) 



where e is the SF efficiency per dynamical time. If e is 
independent of the density, then the result is a Schmidt 
law with exponent ns = 1.5. If, moreover, the scale height 
(H ~ Eg//>g) is also constant, then we have a KS law with 
exponent 1.5, which agrees with the empirically determined 
value at the la level. Note, however, that the assumption 
of a constant scale height is non-trivial. We will show that 
this assumption is unnecessary, generally incorrect, and that 
dropping it gives us considerable insight as well as a much 
more accurate recipe for the implementation of SF laws in 
numerical simulations. 

Numerical simulations o f th e formation of 
galaxies (e.g. ICen fe Ostrikerl Il992b 
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; semi- a n alytic models (e.g. 
Kauffmann, White, fc Guiderdonil Il993h generally use 
a Schmidt law of the form J5}, althou gh other types 
of models have been pr op osed (e.g. iLi et all [2005; 
iBooth. Theuns. & Okamota l2QQ7h and some models 
use different ex ponents for the Schmidt law (e.g. ns = 1; 
lKravtsovll2003h . The Schmidt law is usually combined with 
a threshold volume density, p c , and an upper limit on 
the allowed temperature for SF. Other requirements, such 
as converging flows and Jeans instability are sometimes 
invoked, but these are often sensitive to resolution and tend 
to be relatively unimportant if the problem is numerically 
resolved. The efficiency per dynamical time and the thresh- 
old volume density are the main free parameters, which 
are tuned to make the simulations fit some observable, 
e.g. the observed SF threshold and KS law for the case of 
simulations of isolated galaxies. 

In the next section we will derive the relation between 
the KS and Schmidt laws analytically. We will show that the 
success of the simulations depends on parameters that are 
generally ignored, such as the effective equation of state of 
the multiphase gas and the gas fraction. In section [3] we will 
provide a recipe to reproduce arbitrary KS laws for arbitrary 
equations of state and gas fractions. Contrary to other pre- 
scriptions, the parameters of our module cannot be tuned 
because their values are derived analytically from the ob- 
served surface density threshold and KS law. In section 2] 
we will then use high-resolution hydrodynamical simulations 
to show that our recipe works extremely well. Despite the 
absence of tunable parameters, our prescription is able to re- 
produce observed SF laws with much higher precision than 
has been possible in the past. 



2 RELATING SURFACE DENSITY, VOLUME 
DENSITY, AND PRESSURE LAWS 

The KS law involves surface densities, whereas the Schmidt 
law deals with volume densities. Prescriptions for SF in 
three-dimensional simulations of galaxy formation cannot be 
expressed in terms of surface densities, but do aim to match 
the observed KS law. Although related, there are important 
differences between volume and surface density laws. In the 
literature it is often assumed that the exponents of the cor- 
responding KS and Schmidt power laws are the same, but, 
as we shall see, this is generally not true. 

The KS law describes the SFR when averaged over 
scales that are large compared with individual star clusters. 
In our study of its relation with the Schmidt law, we will 
therefore assume that the volume densities, as well as the 
pressure, are averaged over similar scales. This is precisely 
what is relevant for simulations of galaxy formation which 
currently lack both the resolution and the physical ingredi- 
ents that are necessary to model the multiphase ISM. 



2.1 Thresholds 

Volume and s urface densities can be related as follows 
(|Schavell200ll S04). If self-gravity is important, the density 
will typically fluctuate on the local Jeans scale. For the case 
of self- gravitating discs, this implies that the scale height will 
be of the order the local Jeans scale. Hence, the gas column 
density is of the order of the "Jeans column density" 



(§) 1/2 (/g^tot) 1/2 , 



(6) 
(7) 

(8) 



where 7 is the ratio of specific heats, X the hydrogen mass 
fraction, / = / g //th, with f g the mass fraction in gas (within 
a scale height of the gas) and f t h the fraction of the mid- 
plane pressure that is thermal, and Ptot is the total midplane 
pressure (i.e., including both thermal and non-thermal com- 
ponents) . Putting in numbers yields 



29 M, 



00 i\ /r -2 rl/2 I Ptot k 

28M pc fj ^ 104cm - 3K 



1/2 



(10) 



where we assumed 7 = 5/3, \i = 1.23 and X = 0.752. 

If the gas is far from local hydrostatic equilibrium, then 
the scale height differs substantially from the Jeans length. 
Far out of equilibrium means in this case that either td yn <C 
t sc or tdyn ^> t sc , where t sc = H/c s is the local sound crossing 
time, with H the scale height and c s , e ff = (7-Ptot/ Pg) 1 ^ 2 
the effective sound speed. If left undisturbed, the gas will 
return to local equilibrium on the time-scale min(td yn , tsc), 
which is very short for the densities of interest here (tdyn 00 

1/2. 



1 / VGp « 8 x 10 7 yr / g 1/2 (n H /l cm" 3 ) 1 ) . Note that if the 
entire disc were far from hydrostatic equilibrium, then the 
scale height would fluctuate strongly everywhere, contrary 
to what is observed. 

Following S04, equations (0 and JH} can be used to 
convert the observed surface density threshold for SF into a 
threshold volume density or a threshold pressure. Assuming 
that / - / g - 1 and T - 10 4 K at the threshold for SF, 
and using his prediction for the surface density threshold, 
S c - 3-10 M pc~ 2 , we obtain n H ,c ~ 10" 2 -10 _1 cm -3 
or P to t,c/k ~ 10 2 -10 3 cm~ 3 K. This agrees well with the 
recipes used in numerical simulations, which typically set 
nu,c ~ 10 _1 cm -3 , although some studies of individual 
gal axies have used much higher valu es (e.g., tih,c = 50 cm -3 
for lKravtsovll2003l and 10 3 cm" 3 for lLi et alJl2005h . 

2.2 KS laws 

Using J5J) we can write the gas consumption time-scale © 



(l-n)/2 



(11) 



tg = A" 1 (l M 0P c- 2 ) n (^/ g P tot ) 
Hence, we can write the Schmidt law as 
P* = ^ = A (l Mq pc — 2 ) n (J/ g P tot ) p g . (12) 
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We could write this expression somewhat differently by mak- 
ing use of the ideal gas law, eliminating either p g or Ptot, 
but at the cost of introducing an explicit dependence on T. 
However, the Schmidt law can be further simplified if we as- 
sume that the effective equation of state of the multiphase 
ISM, when averaged over large scales, is polytropic 



Pg 



(13) 



where 7 e ff is the polytropic index, not to be confused with 
the ratio of specific heats 7. This assumption allows us to 
eliminate one variable and we can write, 



p* = A (l M pc 



(n-l) 7eff 

x,lf e P c V"- 1)/2 p g ,e(^ ' 



(14) 



We can see from equation (|14[) that the power- law index 
of the Schmidt law is in general not the same as the power- 
law index of the corresponding KS law. For a polytropic 
equation of state we have 



l)7eff 



n s 



or, equivalent ly, 



+ 1 



2(n s - 1) 

7eff 



+ 1. 



(15) 



(16) 



Hen ce, the two pow er- law indices are equal if 7 e ff = 2 (see 
also ISpringell l2QQ0h , but for any other equation of state 
equality requires n — 1 — ns- 

Simulations of galaxy formation that do not explicitly 
impose an equation of state for gas with a density that ex- 
ceeds the SF threshold, typically predict that the star form- 
ing gas is isothermal, T(p g > p c ) ~ 10 4 K, because the cool- 
ing function cuts off sharply at 10 4 K if molecular cooling is 
not included. Using — 1 gives n — 2 for the commonly 
used Schmidt law index ns = 1.5. To obtain the observed 
value of n = 1.4, one would need to use ns = 1.2. To obtain 
n = 1.4 using ns = 1.5, the simulations would need to have 
7eff = 2.5. 

Equation (fT4)) shows that the Schmidt law depends also 
on the gas fraction, albeit weakly: p* oc /g™ -1 ^ 2 . For n = 1.4 
this gives p* oc / g -2 . If the gas fraction depends on the gas 
density, then this will change the relation between n and ns 
somewhat. If the gas fraction were to decrease with increas- 
ing density, then the effective value of ns would be some- 
what smaller (greater) for fixed n > 1 (n < 1). Conversely, a 
given Schmidt law would correspond to a somewhat steeper 
(shallower) KS law for ns > 1 (ns < 1). 



3 APPLICATION TO SIMULATIONS 

The observed SF threshold, E c , can be converted into a 
threshold for nu or P to t using equations © and J5J), re- 
spectively. If the local gas fraction is not known a priori, as 
will generally be the case for ab initio simulations, then it is 
necessary to assume a constant value. Similarly, the fraction 
of the pressure that is thermal will generally also need to 
be fixed. Following S04, we advise assuming f g — f t h — 1. 
At the threshold, in the outer parts of galaxies, the disc is 



mostly gaseous, there is no SF to inject turbulent energy, 
and the gas is kept warm by UV (background) radiation. In 
fact, the thermal velocity dispersion can by itself account 
for the observed line widths in the outer disc, leaving little 
room for turbulence (S04). When using a volume density 
threshold, we therefore advise using T = 10 4 K, appropri- 
ate for self- gravitating disc with £ g ~ S c (S04). We expect 
the gas to be mostly in atomic form at the threshold and 
therefore advise using the corresponding value of the mean 
particle weight p. 

Only gas with a temperature below some critical value 
T c should be allowed to form stars. For simulations that 
lack either the physics (e.g. molecular cooling and radia- 
tive transfer) or the resolution to model the formation of a 
cold gas phase, critical temperatures in the range 10 4 K <C 
T c < 10 5 K are appropriate. An upper limit on the temper- 
ature is necessary to prevent hot gas from forming stars be- 
fore it cools. This is particularly important for models that 
employ a pressure threshold, rather than a volume density 
threshold, because hot gas can exceed such a threshold even 
for p g <C p g ,c(T = 10 4 K). The results are insensitive to 
the exact value of T c , as long as it is much greater than the 
value below which the cooling time rises sharply (~ 10 4 K 
for most simulations of galaxy formation) and as long as it is 
not much greater than the value for which the cooling time 
is minimum (~ 10 5 K). We use T c = 10 5 K. 

Cosmological simulations should also employ a thresh- 
old density contrast to avoid spurious SF at very high red- 
shift (z> 10 2 ), when the mean density of the universe is sim- 
ilar to the proper density threshold. It is common to choose 
a value similar to the density contrast at the virial radius 
of a collapsed object, A c = pg/p g ~ 60. The precise value 
is unimportant, as long as 10 <C A c < p g , c / pg(z c ), where z c 
is the redshift at which the first resolved objects collapse in 
the simulation. 

The observed KS law, with its normalisation A and 
slope n, can be implemented as the corresponding Schmidt 
law, equation (|12[) , or, if a polytropic equation of state is 
assumed, equation (|14[) . However, since simulations require 
m* rather than /?*, the SFR can be expressed as a function of 
the total pressure, without any direct dependence on either 
density or temperature: 



777,* 



TTTg 



A(lM pc- 2 ) n mg g /gPtot ^ 



(n-l)/2 



(17) 



(18) 



where m g is the gas mass of the element^ for which we want 
to compute m*, and we used (|11[) . This equation can be 
converted to a function of p g , T and p using the ideal gas 
law, or it can be rewritten as a function of p g if a monotonic 
equation of state is assumed. However, by implementing the 
SF law in the form of equation (|18[), the same expression 
can be used regardless of whether an effective equation of 
state is assumed. 

Putting in numbers, we obtain 



1.67 x 10 9 yr 



f S Ptot/k 

10 3 cm- 3 K 



(19) 



2 A gas element can either be a particle or a cell, depending on 
the type of simulation. 
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and 



m* = 5.99xl0~ lu M yr~ 



1 Mr: 



10 3 cm- 3 K 



(20) 



for the observed KS law (equation [2]) and 7 = 5/3. As 
for the threshold, it will generally be necessary to assume a 
constant value of / (note that / g Ptot — fPth)- Fortunately, 
the dependence on / is very weak for the observed value of 
n (ra* oc f' 2 for n = 1.4). 

Simulations must interpret the SF law stochastically. 
That is, the probability that a gas element with SFR p* is 
converted into a star particle in a time step At i^f] 



At 

Prob — min ( — , 1 

tor 



m*At ^ 
mm I , 1 



(21) 



where the gas consumption time-scale and the SFR are given 
by equations JTTJ and ([T5)>. respectively. Since the masses of 
star particles typically exceed the mass of individual stars 
by many orders of magnitude, they are to be interpreted as 
simple stellar populations. Simulations that allow gas ele- 
ments to spawn multiple star particles, rather than convert- 
ing entire gas elements into star particles, should scale the 
probability for SF accordingly: P — > P' = Pm g /m*. 

As an aside, we note that spawning multiple generations 
of star particles per gas element, an approach frequently 
used in the literature, can have undesirable consequences. 
Suppose, for example, that star particles inject energy from 
core collapse supernovae into the surrounding gas particles 
or, in the case of mesh codes, into the nearest cell(s). This 
is a strategy that is commonly used to mimic the effect 
of feedback from SF. As ra*/ra g is decreased, the change 
in the energy per unit mass of the receiving gas elements, 
and therefore their temperature increase, will be correspond- 
ingly reduced. Because the cooling time is sensitive to the 
temperature, this will change the efficiency of the feedback. 
Adding the energy in kinetic form does not alleviate this 
problem: decreasing ra*/ra g will lead to lower velocities and 
thus smaller post-shock temperatures. Kicking the neighbor- 
ing gas elements with a fixed velocity, but with a probability 
proportional to the stellar mass formed also does not solve 
the problem. As ra*/ra g is decreased, the creation of star 
particles would be less and less likely to be accompanied by 
a kick, allowing most of the SF events to proceed without 
any feedback to accompany them. 



3.1 The effective equation of state 

Simulations that do not attempt to simulate star (cluster) 
formation from first principles, i.e., simulations that lack 
the resolution and/or the physics to model a multiphase 
ISM, will benefit from the SF recipe presented here. How- 
ever, since star forming gas (i.e. gas with Ptot > Pc and 
T < T c ) is predicted to be multiphase, one cannot interpret 
the temperature predicted by the simulation for this gas in 



6 Alternatively, one can follow iKatd dl992h and write Prob 
1 — e _At / t s, which cannot exceed unity and reduces to (|21[) for 
At <C t g . What form is preferable depends on the interpretation 
of individual resolution elements. In our simulations the time step 
is always very small compared to t g . Clearly, any simulation for 
which At > t g will produce spurious results. 



the usual way. It therefore also does not make sense to com- 
pute cooling rates, since those depend on the actual kinetic 
temperature of the gas. Instead, it is preferable to specify 
an effective equation of state a priori. Assuming a specific 
equation of state makes explicit what can and what cannot 
be simulated, does not waste resources on calculations that 
are wrong, and also allows one to prevent some numerical 
problems, as we shall see below. 

Although it is possible to devise semi-analytic sub-grid 
models that predict the effective pressure of the multiphase 
ISM as a function of the average densi ty of the gas rep- 
resented by a resolution element (e.g. lYepes etaL 1 Il997l : 
ISpringel &; Hernqu ist 2003), we prefer, for simplicity, to use 
a polytropic equation of state, equation ([T3]) . 

The choice of effective polytropic index 7 e ff, not to be 
confused with the ratio of specific heats 7, can be important. 
Larger values result in a steeper increase of the pressure with 
the gas (volume) density and thus, for n > 1, a shorter gas 
consumption time-scale at a fixed density. This means that 
the maximum densities in the simulation will become smaller 
and the stellar distribution less concentrated. 

For a polytropic equation of state P to t oc pj eii , the Jeans 
length and mass scale as 



1/2 ( 7e ff-2)/2 

4 



Mj OC /g 3/2,(3 7eff -4)/2 



(22) 

(23) 



Hence, 7 e ff = 4/3 gives a Jeans mass independent of the gas 
density and a Jeans length that scales as Lj oc pg 1//3 . If we 
take 7eff = 2, then the Jeans length will be constant, whereas 
Mj oc p g . Hence, for 7 e ff > 2 we do not expect self-gravity to 
promote further collapse, for 4/3 < 7 e ff < 2 we expect gas 
clouds to collapse without fragmenting, and for 7 e ff < 4/3 
we expect both collapse and, possibly, fragmentation. 

A Jeans mass that does not decrease with increasing 
density is desirable from a numerical point of view. The 
alternative, a Jeans mass that decreases with density, will 
lead to a situation in which particles have masses greater 
than the local J eans mass, which is kno wn to lead to artificial 
fragmentation (jBate fe Burkertl 1 19971 ) . A polytropic index 
7eff =4/3 allows one to prevent the strong artifacts that 
result from the inability to resolve the Jeans mass, while at 
the same time allowing collapse to proceed. 

For Smoothed Particle Hydrodynamics (SPH) simula- 
tions 7eff =4/3 not only fixes the ratio of the Jeans mass 
and the mass resolution, it does the same for the ratio of 
the Jeans length and the spatial resolution. Provided the 
gravitationally softening is chosen sufficiently small, the spa- 
tial resolution scales with the SPH kernel h and thus as 
h oc (rag/pg) 1 / 3 oc (ra g /Mj) 1//3 Lj, where ra g is the particle 
mass. Hence, if m g /Mj is independent of the density, then 
so is h/Lj. For =4/3 and a suitable choice of the parti- 
cle mass, the SPH resolution will thus always be sufficient, 
leaving the gravitational softening as the limiting factor. We 
therefore use j e s — 4/3 as our default value, but will also 
explore the effect of varying the polytropic index. 

Note that although our prescription for SF is designed 
to result in a fixed KS law, regardless of the equation of 
state, this does not mean that simulation predictions are 
independent of the equation of state. The effective equation 
of state affects the gas distribution and hence the SF history 
of ab initio simulations of the formation of galaxies. 
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Table 1. Simulations parameters: initial disc gas fraction, / g ; total number of particles, A/tot5 total number of gas particles in the disc, 
^disc5 mass of baryonic particles, m^; mass of dark matter particles, ttzdm; gravitational softening of baryonic particles, e^] gravitational 
softening of dark matter particles, 6dm! effective polytropic index, 7 e fr; KS law exponent, n; gas surface density threshold for SF, S c ; 
initial, central gas surface density, Eq; wind feedback included, (Wind). Values different from the fiducial ones are shown in bold. 
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4 SIMULATIONS 

We ran a suite of 12 numerical simulations of an isolated 
disc galaxy, varying physical and numerical parameters to 
test the theoretical analysis of the relation between the KS 
and Schmidt laws presented in 32 as well as its implemen- 
tation as a recipe for SF as described in 33 Before showing 
the results, we will describe the simulation setup and the 
different runs we carried out. 



4.1 Code and initial conditions 

We used a modified version of the code Gadget 
([Springel. Yoshida. &; White 2001; Springel 2005), which is 
a iV-body TreePM/SPH code. 

We implemented the prescription for SF described in 
section [3] Unless stated otherwise, our models used a SF 
threshold of 7.3 Mq pc -2 , the observed KS law (eq. [2]), 
a polytropic equation of state with index 7 e ff = 4/3 nor- 
malised to T/fi = 10 4 K/0.59 at the thresholcQ, and a gas 
fraction f g = 0.3 which is the initial disc gas fraction of our 
fiducial model. 

A new module for feedback from star formation^, which 



4 We use /i = 0.59 because in our simulations, which include a 
UV background, the gas is still mostly ionized at the SF threshold 
because self-shielding is ignored. Our adopted equation of state 
implies pressure Ptot/k = 2.3 x 10 3 cm -3 K at the SF threshold. 
The value E c = 7.3 Mq pc -2 corresponds to tih =0.1 cm -3 for 
this equation of state and for our fiducial / g = 0.3. 

5 We use kinetic feedback, similar to that of Spri ngel fe Hernquistl 
(2003), except that in our case wind particles are selected locally 
and are not decoupled hydrodynamically. 



will be described elsewhere (Dalla Vecchia et al., in prepara- 
tion), was also implemented, but except for one of our runs, 
feedback was not included in the simulations presented here. 
Radiative cooling and heating was included using tables for 
hydrogen and helium, assuming ionizati on equilibrium in the 
presence of the lHaardt fc Madaul (l200lh model for the z = 
UV background radiation from quasars and galaxies. The 
cooling tables we re generated using the publicly available 
package cloudy (|Ferlandll2000h . 

The initial conditions were generated with a code kindly 
provided to us by Volker Springel. A detailed description of 
the al gorithm is given in lSpringeh Di Matteo, fc Hernquistl 
(|2005h . Here we provide only a brief summary. 

The model consists of a dark matter halo, a stellar 
bulge, and an exponential disc containing both stars and 
gas. The total mass is 10 12 M© (we use h — 0.73) and 
the circular velocity at the virial radius is 163 k ms -1 . The 
dark m atter halo and the stellar bulge both follow [Hernquistl 
(1990) profiles. The Hernquist p rofile is scaled to match 
the in ner density profile of a iNavarro, Frenk, fc White] 
(1996|) profile with concentration parameter c — 9 and 
a mass within the virial radius equal to the total mass 
of the Hernquist profile. The halo has a dimensionless 
spin parameter A = 0.033. The radial disc scale length is 
2.47 h~ x kpc and was compute d by relating it to t he dis c 
angular momentum following Mo Mao, &; White! (|l998h : 
ISpringeL Di Matteo, &; Hernquistl "? 20051 ). The disc contains 
4 percent of both the total angular momentum and the to- 
tal mass. The bulge contains 1.4 percent of the total mass 
and has a scale length one tenth of that of the disc. The 
bulge has no net rotation, while the dark matter halo and 
disc have the same specific angular momentum. The verti- 
cal distribution of the stellar disc follows the profile of an 
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isothermal sheet with a constant scale height set to 10 per- 
cent of the radial disc scale length (i.e. 0.247 kpc). The 
vertical gas distribution is set up in hydrostatic equilibrium 
using an iterative procedure. 

We set the mass of baryonic particles to 



resol -^V nerb 



(24) 



where Mj(p g)C ) is the Jeans mass at the SF threshold for our 
fiducial model (recall that the Jeans mass is independent of 
the density for 7 e ff = 4/3), iV ng b = 48 is the number of 
particles within the SPH kernel and N reso \ is the factor by 
which the Jeans mass exceeds the kernel mass. This choice 
of mass resolution implies 



1/3 



3iVngbmb 
4vrp g 

( 1 \ 1/3 



2/3-7eff/2 



(25) 
(26) 



where h is the SPH kernel and we made use of equations 
(|22[) and (|24[) and we assumed f g to be constant. For our 
fiducial value of 7 e ff this becomes h/Lj = (l/8iV r esoi) 1//3 , 
independent of the density. We chose A^ reso i = 6 which guar- 
antees that we resolve the Jeans mass and (in terms of 
the SPH smoothing) the Jeans length thro ughout the disc, 
easily satisfying the resolution criterion of iBate fc BurkeTtl 
dl997h . Thus, our baryonic particles have mass mb = 5.1 x 
10 4 h~ x Mq. The dark matter particle mass was taken to be 
higher by a factor (Q m — Qb)/^b ~ 4.6. For most of our runs 
the total number of particles in the simulation is 5,000,494, 
of which 235,294 are gas particles in the disc. 

The gravitational softening length was set to 6h = 
10 /i -1 pc for the baryons and to (mdm/^b)^ 3 eb ~ 
17 pc for the dark matter. This is sufficiently small to 
resolve the Jeans length by at least two softening lengths 
up to gas surface densities of ~ 10 4 ' 5-5 Mq pc -2 , which 
greatly exceeds our central surface density of So ~ 2.3 x 
10 2 Mq pc -2 . Note, however, that because our gravitational 
softening can be smaller than the SPH kernel, spurious frag- 
mentation is possible in regions wher e the local Jeans lengt h 
is smaller than the SPH kernel (e.g. [Bate & Burkertlll997h , 
as may happen for polytropic indices smaller than 4/3. 



4.2 Runs 

Table [1] lists all the simulations we have carried out and gives 
the values of all the parameters that are varied. Our fidu- 
cial simulation is labelled fid. Models flO and f90 use initial 
disc gas fractions of 10 and 90 percent, respectively, com- 
pared with 30 percent for our fiducial model. The number 
of particles used in these runs is scaled accordingly, so that 
the particle mass remains constant. Models gammal and 
gamma5 / 3 use effective polytropic indices 7 e ff = 1 and 5/3 
(i.e. isothermal and adiabatic), respectively, rather than our 
default value of 4/3. Run nl.7 uses n— 1.7 as the power- law 
index of the KS law instead of our default value of 1.4. Mod- 
els sigma2.3 and sigma23 have threshold volume densities 
that are an order of magnitude lower and higher, respec- 
tively, than our default value of nu,c = 0.1 cm -3 . Because 
we keep T/fi at the threshold fixed, this yields critical sur- 
face densities that differ by y/l0 from the fiducial value (see 
eq.0). Model wind includes galactic winds, whereas the fidu- 
cial model does not. Our recipe for generating galactic winds 
will be described elsewhere, but we note here that it re- 
sults in much stronge r perturbations to the galaxy than the 
widely used model of ISpringel fc Hernquistl (2003). Finally, 



runs lowres8, lowres64, and lowres512 use particle masses 
that are greater than those used in our fiducial model by 
factors of 8, 64 and 512, respectively. The spatial resolu- 
tion (i.e. the gravitational softening length) is increased by 
factors of 2, 4 and 8, respectively. 



4.3 Results 

4-3.1 Global properties 

Figure [1] shows how our fiducial model evolves. From left 
to right, the panels show face-on projections of the disc gas 
distribution at times 25, 125, and 250 Myr after the start of 
the run. Spiral arms develop rapidly, even though the discs 
were initially axisymmetric. By time T — 250 Myr, the spiral 
arms are starting to fragment. 

Figure [2] compares the face-on projections of the disc 
gas surface density for all runs. Because our galaxies have 
a fixed disc mass, a higher gas fraction implies higher gas 
surface densities, which results in sharper spiral arms. In- 
creasing the polytropic index makes the galaxies more dif- 
fuse, which is not surprising given that both the Jeans length 
(]22[) and mass (|23|) are increasing functions of 7 e ff. On the 
other hand, model gammal, which uses an isothermal ef- 
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Figure 2. Face-on projections of the disc gas surface densities for all models listed in Table [T] All images are 20 h~ x kpc on a side, use 
the same logarithmic colour scale, and correspond to time t = 10 2 Myr. The stellar disc, bulge and the halo are not shown. Each panel 
is labelled by the model ID, as used in Table [T] 



fective equation of state, becomes unstable and fragments 
in sm all clumps. Th e sa me behaviour is visible in the mod- 
els of |Li et alJ (|2QQ5h and lSpringel. Pi Matteo. Hernauistl 
(2005). As we will demonstrate below, this simulation repro- 
duces the desired local KS law, despite the strong instability. 
Increasing the slope of the KS law does not have much effect 
on the morphology of the galaxy. Including galactic winds 
makes the spiral arms more diffuse and results in the gen- 
eration of low density bubbles in the outer disc. Increasing 
(decreasing) the SF threshold results in sharper (more dif- 
fuse) spiral arms. Again, this can be understood in terms 



of the Jeans scales, which decrease with the gas density. 
While model lowres8 still looks reasonably smooth, spurious 
fragmentation is visible in runs lowres64 and particularly in 
lowres512. 

We show in Figure [3] a comparison of the SF histories 
integrated over the entire galaxy from t = till t = 400 Myr, 
which is when we stopped the runs. The SFR of the fidu- 
cial simulation declines gently due to gas consumption. The 
same behaviour is seen in most of the other runs, with the 
exception of the models that we showed to be unstable in 
Fig. [2j gammal and lowres512, as well as f90 which be- 
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t [Myr] t [Myr] t [Myr] 

Figure 3. Comparison of the SF histories of the different simulations. The thick solid line in each plot shows the SF history of our 
fiducial model, fid. The other curves are labelled in the legends of the individual panels. 



comes unstable shortly after the snapshot shown in Fig. [2] 
For these unstable models the SFR increases due to frag- 
mentation, although gas consumption eventually forces it to 
decline. 

Models flO and f90 have initial gas fractions that differ 
by a factor of 3 from model fid. Since the total disc mass 
is fixed, this means that the gas surface densities differ by 
the same factor. Neglecting the change in the total gas mass 
above the SF threshold, which is justified since the SFR is 
dominated by dense gas, we would expect the SFRs to differ 
initially by about a factor 3 1-4 . This corresponds to 0.67 dex 
and agrees nearly perfectly with the simulation results. 

The SF history of model sigma2.3, which has a critical 
surface density smaller than the fiducial model by 0.5 dex, is 
nearly indistinguishable from that of the fiducial run. This 
is because surface densities below the fiducial threshold con- 
tribute very little to the total SFR. As expected, model 
sigma23 is initially below the fiducial one because of the 
higher SF threshold. As the gas collapses to higher densi- 
ties, so does the total gas mass above the threshold surface 
density. The SFR therefore increases slightly, temporarily 
surpassing that of model fid because it is forming stars in 
gas at higher surface densities, until gas consumption forces 
it to decline. 

Run gamma5/3, which has a stiffer equation of state, is 
very similar to model fid. Note, however, that ab initio mod- 
els (i.e. simulations which model the formation of galaxies) 
could still show significant differences between the SF his- 
tories of simulations using different equations of state. This 



would, however, not mean that the predicted KS laws are 
different, which we shall show below is not the case. It would 
merely indicate that the equation of state affects the gas dis- 
tribution. Indeed, run gammal, which has a softer equation 
of state, has a much higher SFR because the gas distribution 
changed as a result of fragmentation. But, as we shall show 
below, the local KS law still holds. 

Decreasing the mass (spatial) resolution by a factor 8 
(2) does not result in a significantly different SF history, in- 
dicating that our fiducial simulation has converged. Decreas- 
ing the resolution further does, however, lead to differences 
because the disc becomes unstable. Again, we shall see that 
this does not imply that the KS laws differ. 

Increasing the slope of the KS law naturally leads to 
a higher SFR, although it eventually drops below the fidu- 
cial rate due to the fast gas consumption. Finally, including 
feedback from SF in the form of winds reduces the SFR. It 
takes some time for the difference to develop because the 
initial conditions were the same, i.e., the starting point is 
not affected by the winds. 

4-3.2 Local properties 

We now come to the most relevant results from the simula- 
tions, the SF threshold and the local KS law. Unless noted 
otherwise, we tested the KS law by projecting the gas den- 
sity and SFR along the minor axis, and measuring their 
values in radial, concentric bins containing a fixed gas mass 
(corresponding to 2000 particles for run fid). For compari- 
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2 gas [M pc" 2 ] E gas [M pc" 2 ] E gas [M pc" 2 ] 

Figure 4. KS law for simulations flO (left panel), fid (middle panel) and f90 (right panel) at t = 100 Myr. Filled circles indicate SFR per 
unit area as a function of the gas surface density, both averaged azimuthally in radial bins containing a a fixed gas mass. The dot-dashed 
lines show the input KS law and threshold surface density for the fiducial simulation. The dotted lines show the expected scaling with the 
gas fraction (see the text for details). Data points with E* < 10 -6 Mq yr _1 kpc -2 are set equal to this value for clarity. The agreement 
between the input SF law and the simulation results is clearly excellent. 
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Figure 5. As Fig. 0] but for simulations gamma5/3 (left panel), gammal (middle panel) and gammal with cylindrical binning (right 
panel). The local KS law reproduces the input law everywhere, but for unstable regions of the disc this is only visible if we use local bins 
rather than azimuthal averages. 



son, we explored two other ways of binning the data: radial 
bins spanning an azimuthal angle Acj) <C 2tt and randomly 
distributed cylinders containing a fixed number of particles. 
All three methods give nearly identical results if the disc is 
stable. All plots in this section are for t — 100 Myr and can 
therefore be compared directly to Figure [2] 

The data points in Figure |4] show the SFR per unit 
area as a function of the gas surface density, both averaged 
azimuthally, for simulations flO, fid and f90. In this figure 
and also in the remaining ones, all bins for which S* < 
10 -6 Mq yr _1 kpc -2 have been set equal to this value for 
clarity. The red dot-dashed lines are not fits to the data 
points, but show the input KS law and SF threshold. 

The agreement between the simulation results and the 
input KS law is nearly perfect, both in terms of its nor- 
malisation and slope. In fact, the agreement is better than 
could have been expected, given that the analytic formulas 
derived in fc|2.2l are no more accurate than the arbitrariness 



in the definition of the Jeans length, roughly a factor ypH. 
However, from equations J7J and JHJ), we can see that this 
would only change the normalisation of the Schmidt law by 
a factor 7r^ n_1 ^ 2 , which is 0.1 dex for n = 1.4, and that the 
slope would not be affected at all. Furthermore, S04 showed 
that for our definition of the Jeans length, J7J) agrees with 
the exact solution for an isothermal, gaseous disc to within 
2 percent. 

The SFR clearly declines relatively sharply below a sur- 
face density threshold that is in good agreement with the 
input value. The fact that the agreement is less good for 
the threshold than for the KS law is not unexpected. For 
a fixed volume density, the corresponding surface density is 
proportional to the Jeans length. Therefore, discrepancies of 
the order of a/tt (0.6 dex) are to be expected, but the differ- 
ence is actually somewhat smaller than that, Note that the 
sharpness of the radial cut off in the SFR is sensitive to the 
way in which the simulation results are binned. Because the 
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galaxies are not axisymmetric, azimuthal averaging smooths 
out the threshold, even if it is sharp locally. 

All three simulations use the same values for the thresh- 
old volume density p g)C and for the normalization of the SF 
law. That is, we always used the fiducial value f g = 0.3 for 
the calculation of p g)C from the input surface density thresh- 
old using equation J7J), and also in the implemented KS law, 
equation ([20]) . We neglected the weak dependence of the SF 
law on the gas fraction because we only know the initial disc 
gas fraction. The gas fraction appearing in the equations de- 
pends on both radius (due to the non-negligible contribution 
of the bulge near the center and the dark halo in the outer 
parts) and time (due to gas consumption). 

The dotted lines in the left and right panels of Fig. |4] 
show the expected SF threshold and KS law assuming that 
the gas fraction appearing in the equations is identical to 
the initial disc gas fraction. As for the uncertainty due to 
the arbitrariness in the definition of the Jeans length, the 
dependence on the gas fraction is only significant for the 
critical surface density and even here it is weak. While run 
flO reproduces the expected threshold better if the gas frac- 
tion is taken into account, the opposite is true for run f90. 
It appears that the initial disc gas fraction does change the 
threshold surface density, but the effect is weaker than the 
naive scaling indicated by the dotted lines. 

Changing the effective equation of state of the star form- 
ing gas has the largest impact on the appearance of the gas 
disc (Fig. [2}. Figure [5] shows the adiabatic (gamma5/3, left 
panel) and isothermal (gamma 1, middle panel) runs, which 
can also be compared with the fiducial model shown in the 
middle panel of Fig. [4] The adiabatic and fiducial runs are 
nearly identical and both agree extremely well with the in- 
put SF law. This shows that we have succeeded in imple- 
menting local KS laws in a way that is independent of the 
effective equation of state. Clearly, in this respect our ana- 
lytic theory describes the simulated galaxies very well. 

At first sight, the isothermal run (middle panel) appears 
to overestimate the SFR by factors of a few for high surface 
densities (E g > 20 Mq pc -2 ). The disagreement is limited 
to the region of the disc that is unstable (Fig. [2]). However, 
what we are seeing here is not a break-down of our theory, 
but rather an artifact caused by azimuthal smoothing. In the 
unstable region the departures from axisymmetry become 
very large and, at a fixed radius, the surface density fluctu- 
ates strongly. Because the KS law is supra- linear (n > 1), 
azimuthal averaging will result in the inference of a steeper 
KS law. 

We can remedy the binning problem by avoiding az- 
imuthal smoothing. In the right panel of Fig. [5] we show the 
results for cylindrical binning. We use an equal number of 
bins as in the other panels, but center them on randomly 
chosen gas particles. For each bin we average S* and E g 
over a cylinder that contains a fixed number of particles. 
We chose bins of 100 particles as a compromise between 
the need to suppress noise (smaller for bigger bins) and the 
need for small spatial bins (using ^> 100 particles per bin 
leads to the same problem as we had for azimuthal binning) . 
Clearly, purely local binning, which is more appropriate but 
more noisy than azimuthal binning, reveals that even this 
manifestly unstable disc follows the local KS law exactly as 
predicted. 

Figure [6] demonstrates that the success of our model 
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Figure 6. As Fig. [4] but for models nl.7 (filled circles) and fid 
(open circles). The prescription for SF can match different input 
slopes of the KS law. 
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Figure 7. As Fig. [4] but for models sigma2.3 (filled circles), fid 
(open circles) and sigma23 (filled triangles). The models match 
the desired SF thresholds. 



is not confined to one particular form of the KS law. If we 
change the slope of the input law from n = 1.4 (open data 
points) to n = 1.7 (filled data points), the simulation still 
reproduces the expected relation. As before, the very small 
deviations at high surface densities are caused by the de- 
crease in the gas fraction due to gas consumption. For model 
nl.7 this happens at a smaller surface density because it has 
a higher SFR (for fixed E g ). 

Figure [7] shows the effect of changing the surface den- 
sity threshold. Comparing runs sigma2.3 (filled circles), fid 
(open circles) and sigma23 (filled triangles), it is clear that 
E c scales as predicted. 

It may be argued that our model, which relies on the as- 
sumption that galaxy discs are self- gravitating, would break 
down for galaxies that are perturbed, although the results 
for the strongly unstable models suggest otherwise. It is 
therefore interesting to see how the model fares if we in- 
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troduce a violent form of feedback from SF. As can be seen 
from Figures [2] and [3] our prescription for winds (which we 
will describe elsewhere) clearly has a dramatic effect on both 
the appearance and the SF history of the galaxy. Neverthe- 
less, as is shown in the left panel of Figure [8] the simulation 
including winds (filled circles) still agrees with the input KS 
law with very little scatter. For surface densities around the 
threshold, the predicted SFR falls somewhat below the input 
KS law. The reason is that at these low densities the wind 
energy is dissipated sufficiently slowly that a significant part 
of the gas surface density is made up of ejected gas, which 
resides far above the disc and has volume densities that are 
too low to partake in the SF. In the right-hand panel of 
Figure [8] we again show models fid and wind, but this time 
we integrate only out to a scale height of 1 kpc. While this 
makes no difference for our fiducial model, it brings model 
wind in nearly perfect agreement with the input KS law. 

Finally, we show the results of a convergence test in 
Figure [9] Reducing the mass resolution by factors of 8 
(lowres8, left panel), 64 (lowres64, middle panel) or even 
512 (lowres512, right panel) does not undermine the agree- 
ment between the simulation results and the input SF law. 
This is quite remarkable, given that the two lowest resolu- 
tion runs suffer strongly from spurious fragmentation (see 
Fig. [2}. Note that the slight discrepancy for the lowest res- 
olution model is most likely due to azimuthal binning (we 
cannot verify this because there are not enough particles for 
spatially resolved binning). 



5 DISCUSSION AND CONCLUSIONS 

Star formation (SF) is poorly understood, but observations 
have revealed the existence of useful laws to describe its 
behaviour when averaged over scales that are large compared 
to individual star clusters. SF is suppressed if the local gas 
surface density falls below S c ~ 3-10 Mq pc -2 , while for 
higher densities a Kennicutt- Schmidt (KS) law, i.e. £* oc 
Eg, is a good description of the data. When averaged over 
entire galaxies, n = 1.4 is a good fit and the same law may 
also hold locally, although this is less well established. If the 
local KS law is supra-linear (n > 1), then the global law can 
only be identical if the total SF rate is dominated by a small 
region, e.g. the nucleus. 

In S04 we provided a theoretical explanation for the 
existence of SF thresholds, showing that the surface den- 
sity threshold for the formation of a cold gas phase agrees 
with the observationally determined SF threshold and that 
this transition triggers gravitational instability. S04 also pro- 
vided a method for the implementation of arbitrary surface 
density thresholds in numerical simulations. Here, we have 
extended the work of S04 to KS laws. Although we have 
not offered an explanation for the observed KS law, we have 
provided a theory that relates various SF laws, that sheds 
new light on previous attempts to account for the observed 
KS law, and that enables the implementation of arbitrary 
SF laws into simulations of the formation and evolution of 
galaxies. 

We showed that the KS law is primarily a pressure law 
since, as is well known, S g oc for a self- gravitating disc, 
although the relation depends also on the gas fraction. The 
KS law is related to the Schmidt law, i.e. p* oc pg S , but 



the relation depends on the scale height and its dependence 
on the surface density. In general, the relation between sur- 
face and volume densities depends on the effective equation 
of state of the multiphase ISM and also on both the gas 
fraction and the fraction of the pressure that is thermal: 
Eg oc /g^ 2 ' f^^T 1 / 2 py 2 . For a polytropic effective equa- 
tion of state, i.e. P tot oc pj eii , the power-law indices of the 
Schmidt and KS laws are related via n = 1 + 2(ns — l)/7eff 
if the gas fraction is independent of the density. More gen- 
erally, the various SF laws, including their normalisations, 
can be related via the equations given in section [2] 

Unless the SF laws are linear (i.e. n = 1 = ns), the 
slopes of the Schmidt and KS laws will generally differ and 
the relation between the two will depend on the effective 
equation of state. Theoretical models have often ignored 
this difference and the dependence on the effective equation 
of state. For example, simulations and semi-analytic mod- 
els typically assume a ns = 1.5 Schmidt law motivated by 
the fact that it implies a gas consumption time-scale pro- 
portional to the dynamical time (p g /p* oc pg 1 ^ 2 ) if, and 
this is sometimes forgotten, the SF efficiency per dynami- 
cal time is independent of the density. It is often mentioned 
that ns = 1.5 is also consistent with the observed KS law, 
n = 1.4 ± 0.1. However, for ns = 1.5 to give n = 1.4 we 
require 7 e ff = 2.5 which is extreme and much greater than 
assumed/predicted by any of the simulations. 

As long as simulations of galaxies lack the physics 
and/or resolution to model the multiphase ISM, in particu- 
lar molecular cooling, radiative transfer, and efficient feed- 
back from stellar winds and supernovae, they will require 
prescriptions for SF that can be calibrated to match ob- 
served scaling laws. We argued that such simulations should 
use an effective equation of state above the SF threshold, 
where the ISM is predicted to be multiphase. A polytropic 
index of 7 e ff = 4/3 is particularly attractive from a numer- 
ical point of view since it results in a constant Jeans mass 
(and, for the case of SPH simulations, a constant ratio of 
the SPH kernel to the Jeans length), while the Jeans length 
still decreases with density. This choice therefore minimises 
resolution effects without directly suppressing collapse. 

Because the observed SF laws involve surface densities 
whereas simulations require volume densities or pressures, 
our work is ideally suited for the implementation of arbitrary 
KS laws into simulations of galaxies. In particular, equation 
(|18[) expresses the star formation rate per resolution element 
in terms of its total total pressure and can thus be used 
even if the polytropic index depends on the density or if no 
effective equation of state is imposed. 

Our equations enable us to predict many of the predic- 
tions of such simulations regarding large-scale SF laws, even 
though our theoretical framework does not explain why real 
galaxies follow such laws. The fact that we can predict simu- 
lation predictions once we know the sub-grid SF implemen- 
tation, implies that simulations of galaxies that do not try 
to actually simulate the multiphase ISM, can only provide 
limited insight into the origin of the observed SF laws. 

We discussed how to implement arbitrary KS laws into 
simulations that do not attempt to simulate the multiphase 
ISM and did so ourselves for the smoothed particle hydrody- 
namics code gadget. We used this code to test our analyti- 
cal relations using high-resolution simulations of an isolated 
disc galaxy. We found that the simulations follow the pre- 
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Figure 8. As Fig. [4] but for models wind (filled circles) and fid (open circles). In the right hand panel E* and E g include only gas 
within a scale height of 1 kpc, while in the left panel we integrate to infinity, as was in the other figures. The model including galactic 
winds matches the input KS law very well, although it falls slightly below the fiducial model for densities around the threshold value. 
Excluding ejected gas (right panel) removes even this small discrepancy. 
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Figure 9. As Fig. HI but for models lowres8 (left panel), lowres64 (middle panel) and lowres512 (right panel). The simulations match 
the input SF law, even for very low resolutions. 



dieted relations nearly perfectly. The surface density thresh- 
old and both the normalisation and the slope of the KS law 
are reproduced with hardly any scatter and without any 
tunable parameters. Moreover, this success is not limited to 
the SF laws that happen to describe the real universe, we 
can reproduce arbitrary input laws, again without tuning 
any parameters. 

The agreement between theory and numerical experi- 
ment remains excellent if we include violent feedback from 
SF or use extremely low resolution. The only significant 
deviations that we found, occurred in unstable galaxies in 
which the spiral arms have fragmented. However, even these 
manifestly unstable galaxies turned out to follow the theo- 
retical relations closely. The apparent discrepancy between 
the input and output SF laws turned out to be an artifact 
of azimuthal smoothing, which becomes inappropriate once 
the galaxy becomes fragmented. 

Finally, we would like to stress that to better constrain 
the models, it is extremely important to get better observa- 
tional constraints on the local SF threshold and the local KS 



law for a representative sample of galaxies. Even azimuthal 
smoothing should be avoided if at all possible, as it greatly 
confuses the comparison between theory and observations. 
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